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Abstract 

Consistency of trans-unification RG evolution is used to discuss the domain of definition of 
the NMSO(10)GUT. We compute the 1-loop RGE j3 functions, simplifying generic Martin-Vaughn 
formulae using constraints of gauge invariance and superpotential structure. We also calculate the 2 
loop contributions to the gauge coupling and gaugino mass and indicate how to get full 2 loop results 
for all couplings. Our method overcomes combinatorial barriers that frustrate computer algebra 
based attempts to calculate SO(IO) /? functions involving large irreps. Use of the RGEs identifies 
a perturbative domain Q < Me, where at the scale of emergence Me < Mpianck the NMSGUT, 
with GUT compatible soft supersymmetry breaking terms, and canonical Kinetic terms, emerges 
from the strong UV dynamics associated with the Landau poles in gauge and Yukawa couplings. 
Due to the strength of the RG flows the Landau Polar regime for gauge and Yukawa couplings 
begins at a scale h.E just above Me and the interpretation of the NMSGUT as an effective theory 
cut off by is well motivated. The large number {Nx) of superheavy fields makes it plausible 
that the observed Planck scale arises {Mpi ~ ^NxMgut) from the wave function renormalization 
of the metric due to GUT scale massive fluctuations but is well separated from it. SO(IO) IR flows 
facilitate small gaugino masses and generation of negative Non Universal Higgs masses squared 
needed by fully realistic NMSGUT fits of low energy data. Running the simple canonical theory 
emergent at Mp through Mx down to the electroweak scale enables tests of candidate scenarios 
such as supergravity based NMSGUT with canonical kinetic terms and NMSGUT based dynamical 
Yukawa Unification. 


2 



I. INTRODUCTION 


Renormalization group equations (RGE) are an important mathematical tool to study 
the evolution of the parameters (couplings and masses) of a quantum held theory with 
energy scale. For example the three gauge couplings of the Standard Model (SM) evolve 
with energy and tend to meet roughly around energy 10^® GeV : this was the hrst dynamical 
hint supporting the “Grand Unihcation” vision [1]-|3]. However the SM has a problem in the 
sensitivity of the Higgs mass to quantum effects of superheavy particles which give rise to 
large loop corrections due to their circulation within loops correcting the Higgs propagator. 
This implies a mass correction : Am\ ~ aM'jr. Supersymmetry (Susy) is the best known 
tool to cure this problem. The two loop RGBs of gauge couplings, superpotential parameters 
and soft terms of a generic softly broken supersymmetric theory have long been available |1] . 
In particular these results give the explicit formulas for the MSSM f3 functions which are 
routinely used to study the evolution of MSSM parameters from UV scales into physically 
meaningful quantities that describe physics near the electroweak scale. 

The combination of supersymmetry and RG flows, given the measured top quark mass. 
Electro-Weak couplings close to those that were hnally measured at LEP(now known with 
better than 1 % accuracy) and the measured strong coupling constant as{Mz) = 0.118 leads 
to nearly exact convergence of the three gauge couplings of the MSSM at = 10^®'^ GeV. 
This striking and robust result has remained the most convincing hint of physics beyond 
the standard model for nearly 30 years since it was predicted to be possible by Marciano 
and Senjanovic |5] if the top quark mass was found to be near to 200 GeV and sin^ 9w was 
larger than 0.23 : as was found to be the case after more than a two decades of searches and 
measurements. Apart from the hints from neutrino oscillations this amazing convergence 
has stood as the unique guide post to the nature of extreme ultraviolet physics for long. 

The closeness of the MSSM Unihcation scale to the Planck scale where gravity becomes 
strong has long tantalized theorists. We have advocated that induced gravity is a natural 
partner for Asymptotically strong GUTs PHO] since their scale of Asymptotic Strength and 
UV condensation should function both as a UV cutoff for the perturbative GUT and set the 
scale for its contributions to the strength of gravity. Recent theoretical arguments [6] renew 
the old speculation |7] that the observed Planck mass may receive dominant or signihcant 
contributions ~ y/NMcur (weakening or inducing gravity by raising the effective Planck 
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mass = Mpi + Mgut) if there are a large number(Atx) of heavy particle degrees 
of freedom of mass Mqut- To our mind the most appealing scenario[2l [ID] is the interpre¬ 
tation of the strong coupling scale of the NMSGUT as a physical UV momentum cutoff. 
Simultaneously gravitational variables(metric,vierbein,gravitino) are demoted to the role of 
a background even as they are supplied with Kinetic terms by the effect of matter quantum 
fluctuations. Their effective action and strength are determined by GUT scale wave function 
renormalization of the dummy variables introduced firstlyto implement general covariance. 
Situational boundary conditions relevant to large scale astrophysical and cosmological con¬ 
texts which are, very plausibly, the only ones where gravity is actually relevant will then 
specify the stress densities that source the classical gravitational fields and waves. Such an 
acceptance of the secondary and induced nature of the gravitational field which needs no 
quantization might hnally lay quantum gravity to rest as an irrelevant incubus, at least to 
the satisfaction of those concerned with testable hypotheses, provided it were anchored in a 
interpretation of the Planck cutoff as a physical cutoff arising from the breakdown of GUT 
perturbativity. Induced gravitational kinetic terms have long been postulated [7] on grounds 
of perturbative wave function renormalization of the graviton due to heavy particles. If we 
take the Planck mass as its experimental value (10^®'^ GeV) and A = Mx = 10^®'^ GeV this 
seems to indicate N ~ 10^. Thus it is very interesting to note that in the NMSGUT there are 
640 chiral superfields and 45 Vector superfields of which 626 are light. This large number of 
SO(IO) coupled helds are precisely what make the couplings diverge strongly in the UV. If we 
count each chiral and vector superheld as 4 degrees of freedom we see that number of heavy 
particle degrees of freedom is V ~ 10^'^. This is in the right ball park to justify the claim 
that the NMSGUT corrections to the graviton propagator actually reduce gravity to the 
weakly coupled theory we observe. By this line of reasoning the Planck scale is determined 
by the unihcation scale of the NMSGUT or its flavour unifying generalization(the so called 
‘YUMGUT’[H] which has even more superhelds). In this picture the Landau Pole(s) of the 
NMSGUT signal a physical cutoff for the perturbative GUT at a scale Ax ~ 

GeV, are the scale of UV condensation driven by SO(IO) gauge forces and moreover set the 
observed strength of gravity. Gonversely the observed strength of gravity actually dictates 
the precise value of the UV momentum cutoff to be used when computing GUT quantum 
effects in any renormalization scheme. Thus the relation between differen cutoff schemes 
is presumably deducible. In the Landau Polar region on the other hand the gauge cou- 
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pling is strong and the theory has entered some sort of condensed phase [9l [TO]. Thus the 
range of scales where the gauge symmetry of the unihed gauge group has unsuppressed play 
seems conhned to a narrow range of scales < Q < GeV. The UV flows of 

Asymptotically free GUTs (of which, in our opinion, no fully realistic example as success¬ 
ful as the realistic Asymptotically Strong Susy SO(IO) models [II], [12] really exists) cannot 
further constrain these scales and only seem to offer the picture of a weakly coupled gauge 
theory crushed as an irrelevance by the strength of gravity above In contrast we 

argue that ASGUTs[ni[in| point to simple yet phenomenologically and calculationally viable 
linkage between gravity and Grand Unihcation of non gravitational forces and matter. 

The very Asymptotic strength of NMSGUT RG flows also hints how the weakly coupled 
gauge theory and a weakly coupled gravitational theory can emerge supernatant at large 
length scales upon the condensate of strongly coupled physics at the smallest length scales. 
The IR flows of these theories very rapidly drive the coupling from arbitrarily strong coupling 
to the weak values observed near the Supersymmetric Unihcation scale gio 95 / rsj Q ^ ^ 
(subscripts 5 and 10 refer to SU{5) and S'O(IO) normalizations for the running gauge cou¬ 
pling constant). From this point of view the trans-unihcation hows of the GUT gauge 
and Yukawa couplings that presumably underwrite the convergence of MSSM couphngs(and 
third generation Yukawa couplings |13]) at or near ~ GeV require the existence 

of a regime Q < Me where a perturbative unihed theory actually operates as the proper 
renormalizable ehective theory describing all particle phenomena except gravity. However 
when -as in the case of the theories studied in this paper- the Susy GUT is Asymptotically 
Strong(AS) in gauge and Yukawa couplings, it becomes problematic even to dehne such 
an energy regime. Nevertheless the nature of the RG hows in the trans-unihcation or sub- 
Planckian regime does have a vital bearing on many interesting physical questions such as 
havour violation in Susy theories [H] and the freedom to choose soft Susy breaking param¬ 
eters required by realistic hts beginning from simple and universal Susy breaking scenarios 
such as canonical Supergravity(cSUGRY) type parameters at the upper limit Me where the 
GUT emerges from the strongly coupled UV regime proper. 

The so called New Minimal Susy GUT(NMSGUT) based on SO(10) gauge group and 
the 210 © 126 © 126 © 10 © 120 Higgs system [TTl IT^ [TSTITT] is the simplest and most phe¬ 
nomenologically successful AS GUT in existence. It has repaid thirty years of detailed 
investigation by exhibiting a remarkable hexibility to accommodate emergent phenomena 
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and their associated data in one overarching calculable theoretical framework and resolve 
long ontstanding problems of nnihcation in terms of the qnantnm effects implied by its 
spontaneons symmetry breaking and associated mass spectra. This has resnlted[Ill |T2] in 
a realistic nnihcation model which is compatible with the known data and with distinctive 
predictions for the Snsy spectra one hopes to observe at the LHC and/or its snccessors. 
Thns it is now topical to examine the RG hows of this theory in the snb-Planckian/trans- 
nnihcation regime to see whether they allow consistent dehnition of a pertnrbative GUT 
over an appreciable energy range, even thongh it is clear that above some scale qnite 
close to the scale of emergence(M^;) of the pertnrbative GUT the theory is in the Landan 
Polar regime where most conplings have diverged and the pertnrbative GUT description 
must be abandoned. 

The NMSGUT requires [TTl [12] small gaugino masses, large squark masses and negative 
non universal Higgs mass squared (NUHM) soft parameters to accomplish EW symmetry 
breaking and ht fermion masses. Such parameters require justihcation, in particular for 
simple cSUGRY scenarios : gravity mediation with canonical gauge and scalar kinetic terms. 
Soft Susy breaking parameters in minimal Supergravity (mSUGRY) are typically assumed to 
be generated well above the GUT scale i.e. near the Planck scale Mp = 

GeV. To consider the ehect of renormalization from Planckian scales to GUT scale, when the 
GUT symmetry is unbroken, one needs the explicit form of GUT RGEs. As is well known 
the NMSGUT exhibits a Landau pole in the generic UV running of the gauge couplingpiTU] 
quite close to the perturbative scale of grand nnihcation. In fact the large coefficients in the 
/3 functions of its other couplings imply the Landau Polar regime involves all couplings. Thus 
there can only be a small energy interval Mx < E < Mp during which the NMSO(10)GUT 
RGEs are usable. Due to the strength of the running it can still have important effects 
even over the short energy range available in ASGUTs as compared to the evolution over 
three decades of energy in the havour violation study of SU(5) SUGRY-GUTs[TT]. If the 
nnihcation program is carried out by running down simple and perturbative data initially 
dehned at Mp using hrst the ASGUT RGEs and then the ehective MSSM RGEs(with added 
neutrino Seesaw and other exotic ehective operators) then the rapid weakening of ASGUT 
couplings towards the IR ensures that a the trans-unihcation how remains perturbative and 
the calculation well dehned. On the other hand the UV how of such theories enters the 
Landau Polar region just above Mp implying that we must assume a physical UV cutoh 
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Ae — Me for the whole Grand unihcation scenario. Beyond this energy lies the true cielo 
incognito where all couplings are no longer weak : “Whereof one cannot speak, one must be 
silent.” 

In spite of their relevance the RGBs for the NMSGUT had so far never been presented. 
In principle the application of the generic formulas of Martin and Vaughn is algorithmic and 
straightforward. However Gomputer Algebra programs [H] that aim to calculate the RG 
functions automatically given the Lagrangian cannot, in practice, handle the combinatorial 
complexity in theories with as many helds as the MSGUT or NMSGUT. Using the vertex 
structure of the superpotential and SO(IO) gauge invariance as constraints makes the sums 
over the components of the large irreps (210, 126 ,126 and 120) required by the formulas 
offfl tractable. The form of the RGBs for supersymmetric theories is governed by the su¬ 
persymmetric non-renormalization theorem [19] whereby holomorphic(superpotential) cou¬ 
plings are free of renormalization except that arising from wave-function renormalization. A 
similar simplihcation is observable in the formulas for the soft couplings and masses. Once 
the tricks for computing the one loop anomalous dimensions are mastered the two loop 
anomalous dimensions and thus f3 functions also follow with some additional combinatorics. 
In this paper we present the NMSGUT one loop /3 functions. However for the case of the 
gauge coupling and gaugino mass we also give the two loop results. We have also calculated 
the two loop results for the rest of the hard couplings and soft Susy breaking parameters |20j 
and we indicate how the methods used for the one loop calculation suffice to yield also the 
two loop results. The other explicit two loop formulas and their effect on running will be 
discussed in a sequel. 

With strict assumptions such as canonical kinetic terms and canonical supergravity type 
soft breaking terms the gravitino mass parameter (^ 3 / 2 ) and the universal trilinear scalar 
parameter Aq (~ 1 ^ 1312 ) are the only free parameters since then there are not even any gaug¬ 
ino masses, the common scalar soft hermitian mass is 7713/2 and the soft bilinear (“B” type) 
parameters are determined by Aq, 7773/2 |2T]. Then the soft Susy breaking parameters at GUT 
scale Mx are determined by running down soft parameters of NMSGUT from Me with just 
these two soft parameters as input. Of course in general one may also consider introducing 
more general soft terms, but our idea here is to show the power of the SO(IO) RG flow to 
generate suitable soft terms at Mx even when placed under such strong constraints. The 
NMSGUT SSB and effective theory are explicitly calculable in terms of the fundamental 
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parameters. However in practice the extreme non linearity of the connection between these 
parameters and the low energy data implies that only a random search procedure (for pa¬ 
rameters dehned at Me) combined with RG flows past intervening thresholds down to Mz 
can hnd acceptable hts of the SM data. The degree of conhdence in the completeness of 
the search diminishes exponentially with the increase in number of fundamental parame¬ 
ters. Thus every reduction in the number of free parameters represents signihcant progress 
towards dehning a falsifiable model. The present work may thus be seen as an attempt 
not only to improve the UV consistency but also to enforce a simplihcation of the htting 
problem by using constraints from the consistency of trans-unihcation RG flows. 

In fact we shall see that the SO(IO) RG flow will identify an additional constraint or 
tuning that must be imposed to keep the soft holomorphic scalar bilinear (‘B’) terms for the 
MSSM Higgs pair in the TeV^ region mandated by NMSGUT htsfUl [12] as well as a RG 
flow based scenario whereby the values of the B parameters may naturally be left in this 
region. Various seemingly peculiar aspects of the NMSGUT parameter choices may hnd an 
explanation in terms of the RG hows at high scales. For instance the negative non universal 
Higgs mass squared parameters ^ which are found in NMSGUT are also justihable by 
the RG hows between Me and Mx- Minimal SUGRY predicts that all soft scalar masses 
squared are positive and equal to fh® scale where they are generated. Soft gaugino 

masses will be generated at two loops from the other soft terms but do not arise at one loop 
if set to zero to begin with. This justihes the typical hierarchy we observed in NMSGUT hts 
whereby sfermions are in the 5-50 TeV range and are much heavier than the gauginos of the 
ehective MSSM (which lie in 0.2-3 TeV range: depending on the lower limits imposed by 
hand in the search). Also the NUHM with negative masses are preferred to have controlled 
lepton havor violation in Susy-GUTs[22|. Similarly choice of the SUGRY emergence scale 
below the Planck scale may also allow adjustment of the gaugino mass and other low energy 
parameters. 

In Section 2 we introduce the formulae of |1] and evaluate them in terms of the parameters 
in the superpotential of the NMSGUT. In Section 3 we present examples of running in 
the sub-Planckian domain. We summarize and discuss our results in Section 4. In the 
Appendix we collect the explicit form of the RG f3 functions of the NMSGUT for soft and 
hard parameters. 


II. APPLICATION OF MARTIN-VAUGHN FORMULAE TO THE NMSGUT 


The generic renormalizable Superpotential without singlets is |lj 

W = (1) 

Here are chiral superhelds which contain a complex scalars 0, and a Weyl fermions ipi. 
The generic collective indices k run over both the different SO(IO) irreps of the NMSGUT 
and dimension of those irreps. The generic Lagrangian corresponding to Soft Susy breaking 
terms is given by 

Lsoftsusy = - ^MAA + h.c. (2) 

jfjk 

are the soft supersymmetry breaking trilinear couplings, bij the soft breaking bilinear 
masses,the Hermitian scalar masses and M is the 50(10) gaugino mass parameter. 
The arrays /i®G 6®-^ are all symmetric and we have allowed for SO(IO) invariant uni¬ 

versal gaugino masses only corresponding to canonical diagonal gauge kinetic term functions 
and SO(IO) invariant 2-loop generation of gaugino masses. 

The theory we now call the New Minimal Supersymmetric SO(IO) GUT was proposed |16j 
by Mohapatra and one of us (GSA) in the early days of supersymmetric GUTs and was 
essentially the first complete and consistent supersymmetric SO(IO) GUT. Its natural and 
minimal structure led another group ini to independently propose it around the same 
time. Following neutrino mass discovery and formulation [23] of high scale left-right and B-L 
breaking Minimal Left Right supersymmetric models (MSLRMs) in the last years of the 
last millenium, it was realized PI in 2003 that it - and not another R-parity preserving 
supersymmetric SO(IO) GUT based on 45 © 54 [25], nor any other competing model such 
as supersymmetric SU(5) with right handed neutrinos added - was the parameter counting 
minimal Susy SO(IO) model. In the same paper it was shown that the GUT SSB can be 
reduced to the solution of a simple cubic equation for one of the vevs. Thereafter it was the 
subject of intense study which calculated its spectra |2SH2S] and specified the roles(couphng 
magnitudes) required of the different Higgs representations for complete fermion fits [301432] 
taking proper account of the role of threshold effects at Ms and Mx(on gauge unihcation). 
Recently we established [121 [IS| that if proper account was taken of the threshold effects at 
Mx on the relation between effective MSSM and GUT Yukawa couplings then the latter 
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-which determine fermion masses and proton decay- can emerge so small as to suppress the 
long standing problematic fast proton decay due to dimension 5 operators completely. The 


210, 126 and 126 Higgs break Susy SO(IO) to MSSM. The 10 and 120 Higgs are mainly 


responsible for the larger charged fermion masses while the small Yukawa couplings of 126 
produce adequately large left handed neutrino masses via the Type I seesaw mechanism, 
instead of failing to do so due to too large right handed neutrino Majorana masses : as feared 
from the early days of this model[T^. Moreover these quantum corrected effective Yukawas 
restore a welcome freedom from the onerous constraints (such as yi, — yr — Vs — Vii, yb/Hr — 
ys/y^J. [S3]) on fermion Yukawas imposed by the NMSGUT proposaljUJ ET] to use mainly 
the 10,120 irreps for charged fermion masses. Thus the model as it stands[IIl [I2] is fully 
realistic and invites further scrutiny. The current study is a part of an effort to simplify 
and unify the fitting procedure by comparing the MSSM parameters implied by randomly 
chosen GUT parameters at to the ElectroWeak and fermion mass data at Mz after RG 
evolution through, and threshold corrections at, the intervening scales. 

The Superpotential of the NMSGUT is ; 

kU = -ynHj + ^^IJKL^IJKL + -^^^IJKL^KLMN^MNIJ + ^^IJKLM^IJKLM 

Tj _ 1 __ 

+ ^^^IJKL^IJMNo'^KLMNO + JKLm{7'^IJKLM + I^IJKLm) 

+ 2 ^;^!) ^IJK^IJK + + ^^Qijk^mnk^ijmn 


1 


2(3!) 


QiJK^KLMn{C^LMNIJ + O^LMNIj) + 




(3) 


Here middle roman capitals /, J, K.... are indices of the vector of SO(IO). All SO(IO) tensors 
are completely antisymmetric and the 5 index ones also obey duality conditions which halve 
their independent components. The indices i, j, k of the generic notation refer to both the 
representation and its internal (independent) components : i = = l....d{R). So 

for example for the 10-plet i = {10;/} , but for the 45-plet i = (45; [U]} with only one 
ordering of each anti-symmetrized pair(J < J) included. Similarly for the 120-plet the 
index i will run over all the 120 different combinations of 3 anti-symmetrized vector indices 
: i = [120; [UK] : I < J < K]. 

As familiar from the MSSM the chiral gauge invariants in the superpotential are the 
templates for the S'O(IO) invariant soft supersymmetry breaking terms. So corresponding 
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to each term in the superpotential we have a soft term in CsoftSusy For example we have 
A corresponding to A, 6 $ corresponding to and a Hermitian mass squared parameter for 
each Higgs representation. In all we have {A, fc, p, 7 , 7 , C) C h, f, g }, &e} and 

{ml, m|, m|, m|), m|-,m|} parameters in the NMSGUT soft Lagrangian,where m| is a 3 
by 3 hermitian matrix. 

Our successful hts mm show that fermion data and EW symmetry breaking requires 
negative Higgs soft masses and soft parameter 6 ^ in the multi-TeV range at GUT 

scale {pH turns positive at low scales) along with gaugino masses in the TeV(gluino) and 
sub-TeV(Bino,Wino) range. In the following sections we will see that such initial values of 
the soft parameters can be generated by running of the SO(IO) theory specihed above even 
over the short range from Me to Mx = Mqut and even when beginning from very restricted 
scenarios for the initial parameter values : such as those implied by cSUGRY. 

We dehne the f3 functions at n-loop order for any parameter x after extracting n powers 
of l/( 167 r^) for convenience in presentation : 

dx _ pi ^ . 

dt (Ibvr^)" 

- The one-loop /^-functions for the SO(IO) gauge coupling and gaugino mass parameter 
M have the generic form: 

71 =- 3C(G)] i /3S>=2/37M/9 (5) 

here S{R) and C{G) are Dynkin index (including contribution of all superhelds) and Gasimir 
invariant respectively. Table [I] gives the Dynkin index and Gasimir invariant for different 
representations of NMSGUT. We get a total index S(R)=l-|-(3x 2)-|-28-|-35-|-35-|-56=161. 
So one-loop P functions for the SO(IO) gauge coupling and gaugino mass parameter are 

= WTsi’o (6) 

7;’ = 2-iMgl„ (7) 

The general form of 1-loop beta function for Yukawa couplings is : 

^ ^ ^ ( 8 ) 

where 7 F) is the one loop anomalous dimension matrix. Thus we need to calculate the 
anomalous dimensions for each superheld. SO(10) gauge invariance implies that 7 * must 
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d 

S(R) C(R) 

= d(G)S(R)/d(R) 

45 

8 

8 

10 

1 

9/2 

16 

2 

45/8 

120 

28 

21/2 

126 

35 

25/2 

126 

35 

25/2 

210 

56 

12 


TABLE I: Dynkin index and Casimir invariant for different representations of NMSGUT 

be field-wise and (irrep) componentwise diagonal. This simplifies their computation enor¬ 
mously. The generic one-loop anomalous dimension parameters are given by 

tL-'= = (9) 

with Yijk = 

To see what is involved in calculating consider the example of the 210-plet. The in¬ 
dependent components of this irrep correspond to non identical combinations of four ordered 
and unequal vector indices : I < J < K < L. Let us select one say 1234. SO(IO) invariance 
requires that 7 / is diagonal so that i = 1234 requires j = 1234 : the propagator correc¬ 
tion will obviously not allow mixing with a different representation than 210. So we are 
required to sum over all possible symmetric combinations of independent 210 components 
pq={{I <J <K < L}, {/' <J'<K'< L'}). 

To calculate we must therefore : 

• Identify the combinations of the chosen component(1234) of <h with other superhelds 
of the model in trilinear gauge invariants. 

• For any given coupling vertex, calculate the number of ways the (conserved) chosen 
(1234) line gets wave function corrections from the helds it couples to in the considered 
vertex. Since it must emerge with the same SO(IO) quantum numbers as it entered 
with and the counting will apply equally to every such held component, a little practice 
suffices to get all 1 -loop anomalous dimensions. 
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Consider first the coupling p^ijkl&ijm&klm 




^4.2.<l>i234 ( 0 12M034Ar “ 013M024M + ©UM© 


23M 


M 

Here M runs over remaining 6 values (M = 5, 6 .. 10 since the 120 plet is totally antisymmet¬ 
ric). In this example we can have 18 possible combinations that couple to ‘hi 234 - Therefore 
the contribution to ( 7 *'^^)i 234 is 


1 

2 


y'f'I’1234-©-©} 


2 


18|pp 

9 


2|pp 


( 10 ) 


Similarly 


7 

-^^^IJKlHm^IJKLM = 74’l234(-f^5Sl2345 + -^6^12346 + 


( 11 ) 


The six allowed index values for H (i.e. 5-10) give-in an obvious shorthand indicating the 
relevant coupling with SO(IO) indices suppressed- 


51 = 6|7p (12) 

H,T, 

The invariant kHjQjKL^ijKL will contribute to 7 ^^^ 
k 

-^Hj&jkL^IJKL = ^4*1234(^^10234 “ ^^20134 + ^^30124 “ H 4 Q 312 ) + .... (13) 

(14) 

H,e 

Thus the anomalous dimension matrix reduces to a common anomalous dimension for each 
independent component of each field and only for the triplicated matter 16-plets need one 
consider mixing. 

In this way one hnds that the one loop anomalous dimension for the 210-plet <h is 
7i'^ = Mk\" + 180|A|2 + 2|p|2 + 240|7p + 6(171' + 171') + 60(1(1' + |CI') - 24^?^ (15) 


Using the anomalous dimensions one can compute the beta functions for all the superpo¬ 
tential parameters. For example the one loop /3 function for A is : 

= 37i.'^A (16) 
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The formulas for the soft terms are closely analogous to those for the Superpotential cou¬ 
plings on which they are modelled. Thus 


+{k i —)■ i) + {k i —>■ i) 

= - 2MY^^^)g‘^C{k) 

+{k i —>• i) -f {k i —)■ j) 

(17) 


where 

7?'^' = 2^1-^nhmnk (18) 

The index patterns of the soft and hard couplings being identical one can calculate one-loop 
/d function for the soft parameter A using the same counting rules used above to sum over 
independent loops. For example the (3 function for the soft trilinear analog of the 210 cubic 
superpotential coupling A (called A) is given by : 

= 3A7?^ + bAyJ^ - 72gl,{X - 2MA) (19) 

where 7 $ ^ and 7 ^^^ = are anomalous dimensions. The hrst was 

given above in eqn(p!5| while its soft (tilde) counterpart is 

= A~kk* + 180AA* + 2pp* + 2A0fjg* + 6 ( 77 * + ^ 7 *) + 60(CC* + CC) ( 20 ) 

These are calculated in the way described earlier with substitution of a soft coupling (h) for 
a hard coupling (Y)(on which h is modelled) and thus the numerical coefficients follow 7 $ 
closely. 

The generic form of the (3 functions for the soft bilinear “B” terms is [1] 

-2{U^ - 2Mp^^)g^C{i) + (i O j) (21) 


Which can again be written in terms of 7 and 7 . Then arguments similar to those given 
above yield : 

13^^^ = 26 ^ 7 ^,'^ + - ASgUh^ - 2Mp^) ( 22 ) 
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Similarly the Hermitian soft masses have generic /? functions 

+ 2Y,„Y>^(mX 

+h^pqyP^ - S5iMM^g^C{i) + (23) 

Again the previous results and a similar one for the doubly soft contribution (i.e. from 
h^^^hipq) yields for example for the 210 soft Hermitian mass : 

/3^^l = 24 '^m| + 720m||Ap + m|(12|7|2 + 12|7|2 + 8|fc|2) 

+m|(8|pp + 120(|CP + ICT) + 8|fcn + m|(480|7p + 12|7p + 120|CP) 
+m|(480|7|2 + 12|7|2 + 120|Cr) + 2^'^ - (24) 

Where 

4 '“ = ( 25 ) 

Thus for example 

7i') = 240|7|2 + 4|/t|2 + 180|A|2 + 2|p|2 + 6(|7|2 + |^|2) + 60(|CP + 1^1') (26) 

As a hnal example of one loop functions consider matter held (Tyi) wave function renormal¬ 
ization due to the matter Higgs superpotential couplings 

w = hAB^Aa(Crj)^^^Bgffl (27) 

where the SO(IO) conjugation matrix C and Gamma matrices Tj may be found in [2H], a, f3 
are Spin(lO) spinor indices and A, B.. are the generation indices. To calculate the contribu¬ 
tion to wavefunction renormalization we need to contract this vertex and its conjugate so 
as to leave as external helds. The remaining numerical factors are : 

hABh*A^B{CTi)a^{CTi)l,p = {C*rjTjC^)^,^{h*h^)A'A (28) 

B,f} 

Thenjin] either C = C[^'^ = Ti x e x ti x e x ti or C = = eXTiXeXTiXe and 

Tj = rj easily give 106a'a{h*h'^)A'A- Similarly the 120plet contributes 120{g*g'^) while the 
126 — 126 pair give 252(/*/^) (since there is a double counting of the 126-plet components 
due to duality within the 252 independent antisymmetric orderings of 5 vector indices). 
Finally since h, f,g are either symmetric or anti-symmetric h*h^ = h^h,g*g'^ = g'^g etc. . 
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A. Two loop anomalous dimensions 


In this paper we study RG flows at one loop level with two important exceptions. Firstly 
the gauge coupling is strongly driven to a Landau pole and it is natural to ask what is the 
two loop correction to the huge positive coefficient in the one loop term. The generic two 
loop formula is 

/Sf = - 6[C'(G)]2 + 2C{G)S{R) + 4S{R)C{R)} - g^Y^^%kC{k)/d{G) (29) 

where the factor in the last term simplifies as G{k)/d{G) = S{k)/d{k). Since for any 
given field type k 'Yhij is diagonal in field type it follows that the sum over k will 

just cancel the dimension of the representation {d{k)) leaving the index S(k) as an overall 
factor weighting the contribution of that field type. This yields 

,3® = 9709gJ„ - + 284” + sG” + 354” + 564” 

+2Tr[4'|) (30) 

The general formula for the two loop gaugino mass (3 function is very similar to the gauge 
beta function 

d” = *9‘{- 61C(9)I" + 2C(,G)S{R) + 4S(R)C(R)}M + - MY“’‘)Yi,t)C(k)ld(G) 

(31) 

and this readily evaluates to 

= 38836^^0^ + 4^?o((7h “ ^4'^) + 2Tr[7^ - + 28(4'^ - + 

35(4'^ - + 35(4^ - + 56(4^ - My^)) 

This concludes the (3 equations we need in this paper. However we have also computed the 
complete two loop results [201 El] • Here we indicate how they are computed. The two loop 
anomalous dimensions y*-^^ are the building blocks of two loop (3 functions and have generic 
form : 

7 ™ = -)r,„„y”'’»K„.F'”7 + 9;„yi„K^«|2C(p)-c(i)| 

+25)94,[C(i)S(fl) + 2C(i)" - 3C(G)C{i)\ (32) 

Again they are field wise and independent component wise diagonal. Only the first term 
requires attention. The intermediate sums over n, r can be broken field wise and thereafter 
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using diagonality of the one loop anomalous dimensions the hrst term collapses to a sum 
over intermediate connected irreps weighted by their one loop 7 s: Thus for example 

Y. Y'^P^V — V I Y. 1 /'QQ'i 

^tmn^ ^pqr^ ■‘■imnn IH ^ ^ imn^ f& ^ J 

Thus the total contribution can be written with the help of one loop anomalous dimension 
parameters. For example : 

7^^^ = -(240|r7|2(4^^ + 7^^)) + 4|fcp(72^ + 7^^^) + + #) 

+360|A|27j,^^ + 4|pp74^^ + 6|7|^(7!y^^ + 7^^^) + 60|C|^(7e^ + 4^^) 
+60|CT(4'^ + 7s ^)) + ^?o(6240|77|2 + 2A\k\^ + 4320|A|2 + 36|pp 
+60|7p + 60|7p + 1320ICP + 1320|Cn + 3864^^0 (34) 

The complete 1-loop anomalous dimensions and (3 functions are given in Appendix. 


III. PROBING THE DEEP CLEFT : APPLICATIONS OF NMSO(10)GUT RG 
EQUATIONS 

A. Landau Polar versus Emergence domain 


Let us begin with the elephants in the room : the huge {3 function coefficients in the 
1-loop gauge /9-function and also in the (3 functions of almost all the chiral multiplet self 
couplings in the superpotential. Thus the coefficients of the cubic terms for the couplings 
{^ 10 , A, 7], 7 , 7 , K, C, C, P} are seen from eqns.(|§ -([ 1 ^ and the Appendix to be (Ibvr^) ^ times 
{137,180,640,142,142,95, 265, 265,16} ! Except for the couplings k, p the other couplings 
grow even faster than the gauge coupling. As noted P, HO] before the huge gauge (3 functions 
imply very rapid change of gio and lead to a Landau pole in the gauge coupling at scales 
within an order of magnitude or so of the perturbative unihcation scale. For the usual (SU(5) 
normalization) value of the gauge coupling at unihcation : = Q;h}(Tf^)/2 = 25.6 

we hnd the SO(IO) gauge coupling has a Landau pole at about A^ ~ ea;p(47r/137a5(M^)) ~ 
10.5M^. In the NMSO(10)GUT, even with the multitude of threshold corrections, a^^(M^) 
can consistently lie in (at most) the range 10 — 40. This corresponds to 9^ = Logio{AE/Mx) 
varying in the range 6^ G [0.4,1.6] although the extreme values are hard to achieve. Thus 
the furthest that one can push the Landau Polar boundary i.e the scale beyond which the 
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theory is certainly fully strong coupled is about GeV. Note however that this ‘UV 

misbehaviour’ pales in comparison with the effect of the combined growth of the Yukawa 
couplings which can reach strong coupling over an scale change by 20% or less ! In fact we 
hnd that this is true of hts found by us earlier when they are extrapolated into the UV. The 
strong divergence and instability of the trans-unihcation flow into the UV implies that it is 
not efficient to look for hts to the complete SM data by parameters thrown at Mx if one 
wishes to have any signihcant range of energies where the SO(IO) GUT is perturbative and 
well dehned. It is quite likely that the parameters optimized for the ht will prove to lead to 
a Landau pole just above Mx- 

Gonversely we propose to turn the strong decrease in couplings in the flow into the IR to 
good account by searching for viable coupling hows valid in the entire range [Mx , Me — A^] 
of dehnition by throwing the parameters -subject to perturbative consistency constraints- 
at a candidate Me and considering downward (i.e. into the IR) running of couplings. The 
scale Me (M^ < Me « Mpianck) is dehned to be the scale where a (weakly coupled) 
ehective GUT with soft Susy breaking has emerged as supernatant to the unknown strongly 
coupled dynamics of trans-emergence scales lying in {Ae, Mpianck)- The strong RG hows, 
as well as the thousands of superheavy particles in the theory with masses ~ Mx make 
a value of Me well below Mpianck plausible without forcing it to coincide with the usual 
unihcation scale Mx- At least prima facie, Me could he anywhere up to GeV. We 

accept that couplings will enter the Landau Polar region at Ae just above Me with the 
reassurance that by choosing the initial values at Me the strongly weakening ehect of how 
to lower scales will reduce the couplings further and ensure that the theory becomes more 
accurately weakly coupled as it approaches the region where the GUT crosses over into 
the low energy ehective theory i.e. the MSSM (with seesaw suppressed neutrino mass and 
other GUT scale suppressed exotic operator supplements). This pattern of energy scales 
is consistent with the expectation [B] that the large number of massive degrees of freedom 
in the NMSGUT will lead to the the Planck mass being dominated by their contribution 
to graviton wave function renormalization : cutoh by the scale of the NMSGUT Landau 
Poles. The Landau Polar latitude Ag is set by examining the UV how from the values 
found to dehne an consistent low energy theory and essentially coincides with Me '- it is 
to be interpreted as a perturbative limit and thus physical cutoh of the ehective SO(IO) 
GUT signalled by the theory itself. It also marks the point where a peculiar and mysterious 
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condensation analogous to confinement in QCD but arising from UV flows takes place in 
the SO(IO) gauge dynamics. In sum, to probe viable scenarios we should take Me to be a 
free along with the hard and soft parameters defined at Me and conduct searches by first 
running these parameters to a matching scale Mx close to the standard MSSM unification 
scale Mx — GeV and then -after applying threshold corrections at that scale |12j 

run the resultant effective MSSM parameters down to the Electro-weak scale Mz and there 
-after applying low scale threshold corrections- match them to the observed standard model 
values. This will be attempted in future improvements of our fitting code. Here however 
we content ourselves by showing that even running couplings down over the short interval 
between the soft Susy parameter emergence scale Me and the GUT matching scale can 
radically reshape the Susy breaking parameter spectrum and bring it closer to the type of 
parameter values we assumed in earlier studiesfUl 1 ^ . 

B. NMSGUT Running down 

To illustrate the actual effect of running down the couplings using the 1-loop jS functions, 
augmented by two loop results for the gauge coupling and gaugino masses, we present an 
example of a flow down from Me = GeV where gw{ME) = 1.0 is quite small enough so 
that RGE flows are perturbative yet large enough that giQ^Mx) — 0.5 required to match the 
MSSM unification value without large threshold corrections can be achieved. We note that 
beginning with ^fio near to 3.0 at the Planck scale(roughly where the RG equation solution 
by 4th order Runge-Kutta methods becomes unstable) one can still flow down close to 
unification gauge couplings ~ 0.5 such as those found for the MSSM unification coupling. 
With lower initial scales the perturbative calculation clearly becomes more palatable and it 
is clear that no purpose is served by starting from a scale which is manifestly unreachable 
given the strength of the flows in the NMSGUT and its cousins. Thus the strong SO(IO) 
gauge RG flows into the infrared, as well as the GUT induced gravity scenario, provide a 
rationale for effective separation of the inevitable gauge-Yukawa condensed strong coupling 
region in the UV from the weakly coupled Susy SO(IO) GUT compatible with the gauge and 
fermion data. As noted before these strong flows can explain and justify certain features of 
the parameter values assumed. 

In our example we take SO(IO) gauge and Yukawa couplings similar to those found in 
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earlier NMSGUT fits using a more rudimentary unification framework where the implications 
of trans-unihcation UV strong RG flows were not taken in to account. Of course only 
an actual search for realistic fits using the strategy proposed here will determine whether 
such couplings can yield a successful fit of the low energy data. We choose soft breaking 
parameters according to canonical kinetic term SUGRY form : all gaugino masses are zero, 
all soft scalar masses equal the gravitino mass at the UV emergence scale: mscaiar{ME) = 
m 3 / 2 (M^) and (for illustration) Aq[ME)=‘ 2^'rnz/2- We also require that the soft bilinears obey 
the strictest form of the gravity mediated scenario| 21 j: 


hi = (Ro - (35) 

at the emergence scale. We chose m 3 / 2 (VfE )=5 TeV and examine the renormalization flow 
from Me to M^. The values of hard and soft parameters at two scales {Me and M^) are 
given in Tables |IT| and III Glearly the RG evolution can be very signihcant and in particular 
the gauge coupling and soft masses change rapidly. Evolution of the Hermitian soft masses 
from Me to Mx is shown in Fig. [T] and we can see that some of them become negative. 
Moreover some of the B parameters also turn negative. In our realistic £ts[lTl [12] we in 
fact find that the values soft hermitian masses squared and B parameter relevant to the 
light MSSM Higgs at the GUT scales need to be negative : which at least the cSUGRY 
framework would contraindicate. 

H, H are constrained to be very light compared to the GUt scale by imposing detT-L = 0 
on their mass matrix (H) which is calculated using the MSGUT vevs[Ill EU 1261129] . The 
left and right null eigenvectors of "H furnish the “Higgs Fractions” UBEllES] whereby the 
composition of the light doublets in terms of 6 pairs of GUT doublets is specified and the rule 
for passing to the effective theory : hi —>■ ajTT, hi —)■ OiiH defined. Then the soft hermitian 
scalar mass terms will give 


mf{h\hi + h\hi) —)■ 


m]j = 


E l |2 2 

|ai| ; 


ml = 


E l - |2 - 2 


(36) 


2 = 1 


Since mj can turn negative when running from Me to Mx we see that negative t^Ijj can be 


achieved. However note that one also has the bij terms for each of the GUT Higgs multiplets 
so that one will in fact also induce the B term for the light Higgs as 


B — hE[(y.\cxi + hE,{(y.2,(y.2 + 0 : 20 : 3 ) + ^$0404 + 60(0505 + ogOg) 


(37) 
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An additional constraint (analogous to that imposed on /x) to maintain the B term at magni¬ 
tudes less than lO^^GeV^ (rather than the RG evolved values which tend to have magnitude 
(Aq ~ 'm- 3 / 2 )l^i >> ''^ 3 / 2 ) ^hus required. The turning of sign of some of the b parameters 
may provide a mechanism whereby the short flow lands these parameters closer to the TeV 
scale values required. 



LogJlOKQ) 


FIG. 1: Evolution of soft masses from Planck scale to GUT scale. Dashed (red), dotted (purple), 
medium dashed (blue), thick dashed (green) and solid (orange) lines represent m^, tuq, m|, 
and m|, respectively. 


IV. DISCUSSION 

We have proposed a framework for a consistent interpretation of Asymptotically Strong 
GUTs(ASGUTs) by considering RG flows of GUT parameters from an emergence scale Me 
of a weakly coupled GUT down to the scale Mx where the GUT is matched to its low energy 
effective theory. Thereafter the MSSM flows from Mx down to Mz determine the experi¬ 
mental predictions of the GUT parameter set chosen at Me- This procedure allows extension 
of the perturbative regime of the unihed theory up to the Landau Polar latitude Ke- Inter¬ 
estingly the large number of degrees of freedom further strengthen the intuition [21 [in] that 
the scale of Gravity may be dominantly set by the effects of (the thousands of) NMSGUT 
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Parameter 

Value at Me = lO^^'^ GeV 

Value at GeV) 

510,55 

1.0, v/2 

0.497,0.703 

A, 5 

-0.0434 -h 0.0078i, -0.313 0.08i 

-0.0133 0.0024i, -0.121 0.031i 

p,K 

0.954- 0.27z, 0.027-FO.lz 

0.21 - 0.06i, 0.0024 0.0088i 

7,7 

0.471,-3.272 

0.0493, -0.425 

c,c 

1.009-h0.831i, 0.36-h0.59i 

0.265-h0.218z, 0.117-F0.192i 

/Ri/10-® 

4.4602 

1.241 

h22/l0-^ 

4.1031 

1.1411 

hss 

0.0244 

.00679 

hu/W-^^ 

0.0 

-1.816-F2.919i 

his/lO-^^ 

0.0 

-1.823-h 1.811Z 

h23/W-^ 

0.0 

-2.955 -h 5.549i 

/ii/io-s 

-.0044-h .16207 

-0.0045 .166z 

/22/IO-® 

6.675 4.8457i 

6.843 -h 4.968i 

/sa/lO-^ 

-9.264 -h 2.7876Z 

-9.498 2.858i 

/12/IO-6 

-0.849 - 1.782i 

-.871 - 1.828i 

/13/IO-6 

.5496-M.1479i 

0.5635-M.177i 

/23/10-^ 

-.4266-F2.231i 

-0.4374 -h 2.287i 

512/10-5 

1.4552 -h 1.599i 

1.016-hi.116i 

513/10-5 

-11.784-F4.9613i 

-8.227-h 3.464Z 

523/10-^ 

-1.6648 - 1.18436i 

-1.162 - 0.827z 

5# 

1015 GeV 

4.55 X lOi^ GeV 

PH 

IO15 GeV 

5.23 X IO 13 GeV 

5s 

IO15 GeV 

5.72 X lOi^ GeV 

50 

1015 GeV 

3.29 X 1014 GeV 


TABLE II: Example of consistent hard NMSGUT-cSUGRY parameters emergent at Me = 

GeV evolved down to GeV using one-loop NMSGUT RGEs for all parameters except 

the gauge coupling and gaugino mass which use two loop evolution. 
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Parameter 

Value at Me = GeV 

Value at M^(10i6-®®) GeV 


-434.0 -h 78.0i, -3127.0 798.0i 

-17.47 + 3.14i, -335.8 + 85.69i 

p,k 

954.4- 269.8i, 273.0-F991i 

-115.7 + 32.7i, -6.39 - 23.19i 

1,^ 

4711 + 0.0i,-32719-FO.Oi 

-80.3 + 0.116i, 142.7 + 0.021i 

c,c 

10091 8305i, 3596 5898i 

125.28 + 103.1Z, 206.77 + 339.14i 

^ 11 / 10-4 

446.02 

63.05 + 0.0028i 

^22, ^33 

4.10,244.19 

0.58 + 2.647 X lO-^i, 34.52 + 0.00158i 


0.0,0.0 

-3.65 + 5.88z, -3.071 + 4.62z 

h23/W-^ 

0.0 

-7.072 + 13.2z 

/ll/10-®,/22 

-0.0436 1.621i, .667 0.4845i 

-0.042 + 1.58z, 0.65 +0.472i 

/33,/i2/10-2 

-9.264 + 2.787i, -0.85 - 1.78i 

-9.013+ 2.71i,-0.83 - 1.73z 

/l3/10“^, f23 

0.55 + 1.15i,-0.427-F2.23Z 

0.535+ 1.12i,-0.415+ 2.17i 

912 

0.146+ 0.16i 

0.073 + 0.08i 

913,923 

-1.178 + 0.496Z, -1.665 - 1.184i 

-0.591 + 0.249i, -0.835 - .594i 

Mg 

0 

-1171.73+ 0.0016i 

6 <j> 

5.0 X IQi^GeV^ 

-3.605 X 10^^ + 6.576 x lO^^^GeV^ 

hH 

5.0 X lO^^GeV^ 

-3.579 X 10^^ + 2.474 x lO^^iGeV^ 

bs 

5.0 X IQi^GeV^ 

3.881 X 10^^ + 6.82 x 10 i 2 ^QgY 2 

be 

5.0 X IQi^GeV^ 

-8.72 X 10^^ - 8.536 x lO^^iGeV^ 

m% 

2.5 X lO^GeV^ 

48070.7GeV2 

mjj 

2.5 X lO^GeV^ 

-1.388 X lO^GeV^ 

ml 

2.5 X lO^GeV^ 

-5.154 X lO^GeV^ 

ml 

2.5 X lO^GeV^ 

1.80955 X lO^GeV^ 

ml 

2.5 X lO^GeV^ 

9.564 X lO^GeV^ 

Eval m| 

>1' 

2.5 X lO^GeV^ 

{2.7892,2.7892,2.7889} x lO^GeV^ 


TABLE III: Values of NMSGUT soft parameters at two different scales evolved by using one-loop 
SO(IO) RGEs. Aq = 10 TeV, m ^/2 = ^ TeV. 23 



superheavy particles. Thus Mpi can he well above Ap which nevertheless plays a part in 
raising Mpi by serving as the physical cutoff scale for graviton wave function renormaliza¬ 
tion corrections due to the NMSGUT as well as the scale for SO(IO) “conhnement” [S]. We 
presented the NMSGUT RG equations to determine the RG evolution of couplings between 
the scale(ME) where the perturbative effective theory (NMSGUT plus weakly coupled and 
softly broken = 1 supergravity) emerges and the matching scale between GUT and the low 
energy effective theory (i.e the MSSM) at Mx- To illustrate the application of these results 
we evaluated the effects of the 1-loop evolution on randomly chosen sets of parameter values 
assuming a minimal, canonical kinetic term, supergravity scenario for the starting parameter 
ansatz. From the Tables and Fig. 1 we see that the RG evolution has dramatic effects on 
the soft Susy breaking parameters. Firstly most of the soft Susy breaking Hermitian masses 
squared of the SO(10) Higgs irreps become negative even though they start from a common 
positive mass. This provides a potentially robust justification of the negative values of ^ 
used in NMSGUT hts [TTl |T2]. Note that the distinctive normal s-hierarchy at low scale is 
strongly correlated with the large negative g we use in the hts. Gaugino masses(MA) 
will be generated by two loop RG evolution between Mx and Mz, even if Ma= 0 at the scale 
Mx- On the other hand the same applies to the evolution between Mp and Mx- Thus even 
canonical gauge kinetic terms in the GUT can still generate adequate gaugino masses. This 
pleasing since we have always resisted invoking non canonical Kahler potential and gauge 
kinetic terms on grounds of minimality/predictivity and to preserve renormalizability of the 
gauge sector. 

Another notable effect is the intermediate scale {0{m^/2Mx)) values of the soft parame¬ 
ters required by the canonical mSUGRY ansatz and induced by the dependence ^ 

~ Mx'ms/ 2 - So we may need to impose an additional condition in order that the contribu¬ 
tion from the soft terms to bp p is 0 ( 7713 ^ 2 ) unless this is achievable via the RG how of Bg 
towards negative values itself. 

The running of trilinear soft coupling and s-fermion mass squared parameters (m|) will 
give distinct values at the GUT scale for the three generations (considered same in earlier 
studies of NMSGUT [TT[ [12]). Our illustration was computed at 1-loop only for simplicity. In 
sequels we will integrate these RG hows with our previous code that incorporates the MSSM 
hows between Mx and Mz- Then one will throw the core soft parameters 7 / 73 / 2 , Aq at Ax 
and run down over thresholds to Mz with one additional hne tuning constraint. Thus the 
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total number of soft parameters will be significantly reduced. Finally the straightforward 
(since the superpotential vertex connectivity is preserved) generalization of these results to 
the case of YUMGUTs[H] will allow us also to perform the RG flows from the Planck scale 
for dynamical flavour generation models based on the MSGUT. These theories have around 
6 times as many helds as the NMSGUT and are thus even more capable of separating Mx 
and Mpi. We note that the techniques we have used to actually evaluate the 2-loop RGBs 
have overcome the combinatorial complexity that prevented their calculation by automated 
means. They can be used for any Susy GUT. 


25 


Appendix 


One-loop RGEs 


One-loop anomalous dimension parameters associated with different superfields: 


- 2gl,5lC{i) 

7^'^ = 200|7|2 + 10|7p + 100|Cp-25(72^ 

= 200|7p + 10|7r + 100|Cr + 32Tr[/t./]-25(7?„ 

4 ^ = 84|«;|2 + 126(|7r + |7r) + 8Tr[ht.h]-9(72„ 

7^'^ = r(|«:p + IpH + 105(|CP + ICT) + 8Tr[^?t.^] - 21gj, 

= 252/t./ + UOglg + lOh^.h - ^ (38) 

1 

_ _Y. wiP? 

H ~ 2 

4'^ = 200|7|2 + 10|7p + 100|Cp 

= 200|7p + 10|7r + 100|CP + 32Tr[/t./] 

= 84|/t;P + 126(|7pIdP)-h 8Tr[/i^/i] 

7^'^ = 7(|«:r + |pp) + 105(|Cr + |CT) + 8Tr[<7t.^] 

= 252/^/ + 120glg + lOh^ .h (39) 

One-loop beta functions for the SO(IO) superpotential parameters and Yukawa couplings 
are: 


/9l'^=37i'^A ; 

= ’K7e ’ + 7e ’ + 7*’) 

(40) 

7(7h^+ 7s^+7#^) ; 

4“’ = 7(7U’ + 7e ' + 7»’) 

(41) 

/c(7h^+ 7e^+7#^) ; 

7'=C(7i'>+77+77) 

(42) 

= C(7e^ + 7g ^ + 7$ ^) 

; tf’=M77 + 27W) 

(43) 
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(44) 


= gi^e - i'yS^f-9 + 9-lS^ (45) 

~ ^7$ V^* ; = ‘^Ih^ f^H (46) 

= (7s ^ + 7g Vs ; Vj = 27©Ve (47) 


Soft parameters RGEs 


;yWj ^ iy. Upq 
h 2 

= 20077* +1077* + loocc* 

= 2OO77* + IO77* + lOoV + 32Tr[/^/] 

= 8 AkK * + 126(77* + V) + 8 Tr [ h ^. h ] 

41 ) = 7(««* + pp*) + 105(CC +V) + 8Tr[(7t.^] 

V = 252/^/ +120^^^ +lO/i^/i (48) 

= lh,pqh^^^ 

V = 240|r)|2 + 4|«|2 + 180|A|2 + 2|p|2 + 6(|7p + |7r) + 60(|Cr + lV) 

V = 200|r)|2 + 10|7r + 100|Cr 

V = 200|7|2 + 10|^|2 + 100|V + 32Tr[/t./] 

V = 84|Kp + 126(|7|2 + |^n + 8Tr[/it./i] 

V = 7(|«p + |pn + 105(|C? + lV) + 8Tr[r.^] 

= 252/^/ +120^^^ +lO/i^/i (49) 

One-loop beta functions for the soft parameters: 
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4^^ = vili^ + + 7g + 2r/(7$ ^ + 7s ^ + 7g ^) - 7%io(7 “ 2Mr/) (50) 

= 7 ( 7 $ ^ + 7s ^ + 7^^) + 27 ( 7 $ ^ + 7s ^ + I^h) - 58 ^ 10(7 - 2 M 7 ) (51) 

= ^( 7 $ ^ + 7g ^ + 72^) + 27 (7^^ + 7^^ + 7^2 - 58^10(^ - 2 M 7 ) (52) 

= «(72^ + 7^^ + 7^2 + 2^(7^^ + 70^ + 7^2 “ 54g^o(^ - ^Mk) (53) 

= ^(7^^ + 27^2 + 2p(72^ + 27^2 “ 66g^o(P “ 2Mp) (54) 

= C(72^ + 7e^ + 70 ) + 2C(72^ + 7^^ + 7^2 - 70pio(C - 2MC) (55) 

= ^( 7 ^^ + 7^^ + 7^2 + 2C(72^ + 7^^ + 70 ) - 70pio(^ - 2MC) (56) 

+ /i. 72^ + (7i^)^-h + 2jjJ^h + + (722^./i) - yP?o(^ - 2M/i) (57) 

7^^^ - ^-72^ + (722^-^ + 270^.g + 2(p.72^ - (722^-^) - y^?o(^ - 2Mp) (58) 

7g 7 + (722^-/ + 27g V + 2(/.72^ + (722^-/) “ y^?o(/ “ 2M/) (59) 

= 26$72^ + 4/i<i.72^ - 48^10(^< 1 - “ 2Mp$) (60) 

= 26^/72^ + 4pB7iJ^ - 18Pio(bii - 2 Mph) (61) 

j3le = + 4fX07^0^ - 42gl(b0 - 2M/ie) (62) 



— ^s(7e ^ + 7g + 2/iE(7s^ + 7g - ^^gioipY, - 2M/is) 


(63) 


= 24')m| + 720m||A|2 + m^(12|7p + 12|7|2 + 8|«;n 

+m|(8|pr + 120(|CP + ICr) + 8|«:r) + m|(480|7p + 12|7r + 120|Cn 
+m|(480|7p + 12|7|2 + 120|Cn + 27^'^ - ^m?9lo (64) 

= 27fJ^m^+ m|(252(|7p + |7p) + 168 |kP) + 168m||fi;p + 252m||7p 

+252m||7|^ + 27^;^ - 36|M|Vo + 32Tr[/i^m|./i] (65) 

/^m| = 27e^m| + m|(14(|K|2 + \p\^) + 210(|CP + ICH) + 14m||pp + Um]j\K\^ 

+210m||Cr + 210m||CP + 24'^ - 84|M|2p2^ + 32Tr[pt.m|.p] (66) 

= 24')m| + m|(400|7|2 + 20|7r + 200|Cr) + 200m||Cp + 20mi|7r 

+400m||7|2 + 24'^ - im\M\^gl^ (67) 

= 24)771!+ m|(400|7|2 + 20|7r + 200|cr) + 200m||cp + 20mi|7r 

+400m!|7p + 24^ “ 100|^p5'io + 128Tr[/'*'.m|./] (68) 

42 = + "^1-4^ + 10/7^m|./7 + 1205(^m|.5( + 2h2f.m\.f + lOmjjh^ .h 

+ 120mlglg + 252mlflf + 24^ - A5\M\^gl (69) 
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